home *** CD-ROM | disk | FTP | other *** search
/ Shareware Grab Bag / Shareware Grab Bag.iso / 008 / 68kgraph.arc / SQRT.ASM < prev   
Encoding:
Assembly Source File  |  1986-12-13  |  1.4 KB  |  39 lines
  1. ******************************************************************
  2. * COPYRIGHT (C) 1986 by Donald Krantz and James Stanley
  3. * - Note: This is a real, live, actual, registered copyright,
  4. *   and should be treated as such. This source code is from
  5. *   the book "68000 Assembly Language", Krantz and Stanley,
  6. *   Addison-Wesley Publishing Company, Reading, MA, 1986.
  7. * 
  8. *   Permission granted by the authors for non-commercial use
  9. *   in programs released to the public domain, as long as this
  10. *   copyright notice remains attached and visible.
  11. *
  12. *****************************************************************
  13. * SQRT - Integer floored square root
  14. *    Entry:     D0.L = number to take root of
  15. *    Exit:    D0.L = square root
  16.  
  17.     xdef    sqrt
  18.  
  19. sqrt:
  20.     movem.l    d1-d3,-(a7)
  21.     clr.l    d1        * D1 will be test case root
  22.     move.w    #$8000,d2    * D2 is succ. aprox. mask
  23. lp1:
  24.     or.w    d2,d1        * load in mask value
  25.     move.w    d1,d3        * d3 is a scratch register
  26.     mulu    d3,d3        * make a try
  27.     cmp.l    d3,d0        * take a look - how'd we do?
  28.     bhi    sk1        * try less than goal - keep bit
  29.     not.w    d2        * we need to mask bit out 
  30.     and.w    d2,d1        * take out trial bit
  31.     not.w    d2        * fix mask
  32. sk1:
  33.     lsr.w    #1,d2        * move mask down
  34.     bcc    lp1        * mask mot done - go around.
  35. sq_exit:
  36.     move.l    d1,d0        * transfer answer
  37.     movem.l    (a7)+,d1-d3    * restore registers
  38.     rts
  39.